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1 Introduction 



The Bhabha scattering, e"^e~ — > e"*"e~, plays a special role in the phenomenology of 
particle physics, since it is employed to determine the luminosity of e^e~ colliders 
operating at both high (~ 100 GeV) and intermediate (~ 1 — 10 GeV) center of mass 
energies. The luminosity of such facilities is essentially determined by calculating 
the (inverse) ratio of the theoretical Bhabha scattering cross section to the number 
of Bhabha scattering events observed. The Bhabha scattering is a process particu- 
larly suited for luminosity monitoring purposes, since its cross section is large and 
QED dominated (this is true in specific kinematic regions: the small scattering angle 
region for what concerns high energy colliders, the large scattering angle configura- 
tion for machines operating in the ~ 1 — 10 GeV center of mass energy range). As 
a consequence, this process can be calculated in perturbation theory with a high 
degree of accuracy. Since the theoretical error on the Bhabha scattering cross sec- 
tion affects directly the precision of the luminosity of an e^e~ collider, a remarkable 
amount of work has been devoted to the study of the radiative corrections to this 
process, both virtual and real, in the last three decades. For an extensive list of 
references on the subject we refer the interested reader to the review in jTj. 

In this paper, we focus our attention on the technical problem of the diagram- 
matic calculation of the virtual corrections to the Bhabha scattering cross section. 
The one-loop radiative corrections to this process, together with the corresponding 
real corrections involving the emission of a soft photon, calculated within the frame- 
work of the full electroweak Standard Model, are well known [2 . At the two- loop 
level, the calculation of the virtual corrections becomes much more involved. Even in 
pure QED, the complete (non approximated) two-loop differential cross section, i. e. 
the two-loop amplitude interfered with the tree level matrix element and summed 
over all spins, has not yet been calculated. A large amount of work was devoted to 
the study of the contributions enhanced by factors of ln(s/mg) [3]. In recent years, 
the complete two-loop QED virtual cross section has been obtained in the limit of 
zero electron mass TJ. The IR divergent structure of this result has been studied in 
detail in 

The main technical challenge encountered in calculating the complete, non ap- 
proximated, set of two-loop virtual correction is represented by the evaluation of 
the box-diagrams. Recently, we proved that the calculation of the two-loop box 
diagrams involving a closed fermion loop can be reduced to the calculation of a set 
of 14 scalar Master Integrals (Mis) 6J. This result has been obtained by means of 
the (by now standard) Laporta algorithm [7j, for the systematic exploitation of the 
Integration-By-Parts (IBPs) |H|, Lorentz invariance (Lis) [HI, and general symmetry 
PU] identities. The Mis involved in the calculation have been evaluated in P CHI , 
within the context of D-dimensional regularization scheme |TTj, and by employ- 
ing the method of the differential equations in the external kinematical variables 
\12\ ITHl E] • The result was conveniently expressed in terms of 1- and 2-dimensional 
harmonic polylogarithms (HPLs, 2dHPLs) ^3 ^21 El CHI CSl ? of maximal weight 3. 

The results in [HI completed the calculation of all the Mis that are necessary 
in order to evaluate the two-loop Bhabha scattering differential cross section in 
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pure QED, if we limit our analysis to diagrams involving a closed fermion loop 
(conventionally referred to as the a'^{Np = 1) differential cross section). In fact, the 
relevant vacuum polarization diagrams have been calculated in [20] • The two- loop 
vertex corrections involving a closed electron loop have been obtained in 122] , 
where the authors present a complete calculation of the vertex form factors in QED, 
without neglecting the electron mass. 

The purpose of the present paper is to provide the expression of the contribu- 
tion of the virtual corrections to the Bhabha scattering differential cross section at 
order a'^lNp = 1) in QED. The UV renormalization of all the two-loop diagrams 
is carried out within the Z)-continuous dimensional regularization in the on-shell 
renormalization scheme. The residual IR divergence in the final result appears as a 
pole in (D — 4), to be canceled, as usual, by the corresponding divergence due to 
the real emission of soft photons in the complete expression of the physical cross 
section. 

The paper is organized as follows. In Section |2] we describe the various steps of 
the calculation, introducing our notations and conventions. In Section|2]we recall the 
main points of the UV on-shell renormalization scheme. In Section0]we re-derive, for 
completeness, the one-loop virtual corrections and the Bhabha scattering differential 
cross section at order a^. The main result of this paper, the virtual differential cross 
section at order a^lNp = 1), is presented in Sectional In Appendix \K\ we collect 
the explicit expressions of the auxiliary functions introduced in the paper in terms 
of HPLs, and we briefly review the analytical continuation of the HPLs from the 
Euclidean to the physical kinematic region. In Appendix El we give the expansions 
of the auxiliary functions in different kinematical regimes. 



In this Section we summarize the strategy followed in carrying out the calculation. 
In order to fix the notation and conventions employed throughout the paper, we 
start by discussing the tree-level Bhabha scattering differential cross section. 
We consider the following scattering process mediated by photon exchange: 



The incoming electron, incoming positron, outgoing electron and outgoing posi- 
tron have momenta pi,P2,P3 and p^, respectively, and finite mass m. The process is 
well described in terms of the Mandelstam kinematic invariants s, t and u, defined 
by the relations: 



2 The Bhabha scattering 
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Figure 1: Tree-level diagrams. 

and where E is the particle energy in the center of mass frame of reference and 9 is 
the scattering angle in the same frame. 

The differential cross section for the Bhabha scattering process can be written, 
in perturbation theory, as an expansion in the fine-structure constant a: 

da{s,t,m'^) (icro(s, t, m^) /a;\ (icri(s, t, m^) /ax"^ d<72{s,t,m'^) ( fa' 

dn " dn dn "^wJ dn ^ 

(6) 

where cro(s,t,m^) is the tree-level (Born) cross section and (Ji{s,t.,m?) {i = 1,2) are 
the higher order contributions. 

The s- and t-channel diagrams contributing to the scattering process at the tree 
level are shown in Fig. ^ where, as usual, the wavy lines represent photons and the 
arrows on the fermionic (solid) lines follow the flow of the negative charge. 

The tree-level amplitude, J^q, is given by 

Mo = At-As, (7) 

where the matrix element of the t-channel (s-channel) diagram in Fig. is indicated 
by At (As). In Eq. (|7j) we indicated explicitly the reciprocal negative sign between 
the s- and t- channel diagrams; in the following, the expressions of the contribution 
of the s-channel virtual diagrams to the differential cross section are given including 
this extra minus sign. 

After averaging the squared amplitude at order over the spins of the initial 
states and summing over the spins of the final states, the differential Born cross 
section reads: 



d(TQ{s,t,m'^) 



dQ s I s 
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st + — + {t- 2m^f 



1 



st + — + {s- 2m^f 



Ob \ 



The first, second, and third term between curly brackets represent the contribution of 
the squared s-channel diagram, the squared t-channel diagram, and the s-t-diagram 
interference, respectively. 

If we consider only the virtual contributions at order a^, it is necessary to evaluate 
the diagrams shown in Fig. |21 The corresponding squared matrix element, summed 
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over the spin of the initial particles and averaged over the spins of the final ones, will 
be given by the interference of each diagram in Fig. |21 with the tree-level diagrams 
of Fig. ^ The differential cross section assumes the following form: 



spin ■■ ' ' ^ — - — 

where the superscript V stands for "virtual" (and the sum is performed on all spins). 
The contributions of the various terms on the r.h.s. of Eq. @ to the differential 
cross section will be discussed in Section HJ 

At order a^, if we limit our analysis to the diagrams containing a closed elec- 
tron loop {Np = 1), it is necessary to consider the contributions coming from the 
interference between the diagrams shown in Fig. 01 and the tree-level amplitude, as 
well as the interference between the one-loop vacuum polarization diagrams with 
the other one-loop diagrams without a fermionic loop. Diagrammatically, we have: 



q;\ 2 daXis, t, m^) 

71 ^ 




+ c.c. 



+ \)-KK X +C.C. + (10) 



The expression of the various contributions to the differential cross section shown 
in Eq. ()10|) will be given in Sectional 
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The evaluation of the virtual radiative corrections to the Bhabha scattering, 
Eqs. fl9ll0|) . leads to the calculation of UV- and IR-divergent integrals, which need 
to be regularized. Both classes of divergences are regularized within the framework 
of Dimensional Regularization scheme fl], in which the regularization parameter 
is the dimension of the space-time D. In order to evaluate the virtual radiative 
corrections we follow a procedure consisting of three main steps: 

• First, we consider the Dirac structure of the squared matrix elements; the 
sum over the spins generates traces over the Dirac indices, that can be easily 
evaluated in D dimensions. 

• At this stage, the squared amplitude describing the interference of a one- 
loop or two- loop diagram with the tree-level matrix element, as well as the 
interference of a one-loop with another one-loop diagram, can be expressed as 
a sum of a large number of scalar integrals (several thousands). According to 
OCni! it is possible to express all these integrals as a combination of 14 Master 
Integrals. This result was obtained by using the Laporta algorithm [7j in order 
to exploit the information contained in various sets of identities (Integration- 
by-Parts, Lorentz invariance and other general symmetry identities) relating 
the scalar integrals among themselves. Once the reduction to the Mis has been 
carried out, we can expand the squared amplitude, which will depend on the 
dimensional regulator around D = A. Using the expressions for the Mis 
also given in [HI Ej, we finally obtain an analytic result, expressed in terms 
of 1- and 2-dimensional harmonic polylogarithms, in which both UV and IR 
divergences are regularized by the same parameter D and appear as poles in 



• The last step concerns the UV renormalization. By subtracting the appropri- 
ate counterterm graphs, the poles in [D — 4) corresponding to the unphysical 
UV divergences of the diagrams cancel out. We perform the subtraction of 
the UV poles in the on-shell renormalization scheme. The result obtained in 
this way still contains IR poles in (Z) — 4), which are canceled in the physical 
cross section by the contribution of the real soft photon radiation. 

3 Renormalization 

We perform the UV renormalization in the on-shell scheme. The counterterm dia- 
grams for the renormalization of the one-loop and two-loop cross section are shown 
in Fig. m 

The counterterms give rise to the following Feynman rules: 



(D-4). 
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Figure 4: Counterterm diagrams at one and two loops. 
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The one- and two-loop renormalization constants are: 



1 + (Z) - 4) + O ((D - 4)2) 



2{D-A) 

zi^^\D) = z[^^\d) (Ward identity) , 



(13) 
(14) 

(15) 
(16) 



(17) 
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(20) 



m 




1 



947 



2C(2) + 0(D-4). 



(21) 



8(D -4) 



+ 



288 
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+ 0{D - 4) . 



(22) 
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The renormalization constants are defined in such a way that the UV-renormahzed 
quantity is obtained by adding to the UV-divergent diagrams the corresponding 
counterterm graphs. 

At the one-loop level, it is necessary to renormalize the photon vacuum polar- 
ization graphs, shown in Fig. |2 (a) and (b), and the 7//- vertex graphs, shown in 
Fig. 121 (c)-(f). In the first case, we add the UV-divergent diagram of Fig. |21 (a) (or 
(b)) and the corresponding counterterm of Fig. |3] (a), ending up with the renor- 
malized expression for the photon self-energy given in the next Section. For the 
one-loop vertex it is necessary to add the UV-divergent diagram of Fig. |21 (c) (or 
(d)-(f)) with the corresponding counterterm of Fig. 0] (b). The UV-renormalized 
form factors found in this way are given in Eqs. ()85I86|) . 

In the two-loop case, we have to renormalize the UV divergences occurring in 
the diagrams of Fig. El 

The three two-loop vacuum polarization diagrams have to be added with the 
corresponding seven counterterm diagrams shown in Fig. 01 (c)-(h) and (1). Note 
that, because of the QED Ward identity, Zi{D) = Z2{D) and the contributions of 
diagrams (c), (d), (g) and (h) in Fig. ^ cancel among themselves. The renormal- 
ization of the two-loop vacuum polarization consists, therefore, only in the one-loop 
electron mass renormalization and in the subsequent cancellation of genuine two- 
loop UV divergence against the counterterm diagram (1). For the two-loop vacuum 
polarization vertices. Figs. El (g)-(j), we have to add the counterterms (i) and (k) 
of Fig. m the first one cancels the sub divergences due to the one-loop vacuum po- 
larization insertion, while the second cancels the overall two-loop divergence. The 
box graphs, shown in Figs. |21 (k)-(r), do not present two-loop UV divergences and 
their subgraph divergence is renormalized by adding the corresponding counterterm 
diagrams of the kind shown in Fig. IH (j). Finally, the diagrams in Fig. |21 (s)-(v) are 
renormalized by adding to them the corresponding counterterm graphs subtract- 
ing the one-loop divergences of the vacuum polarization insertion and of the vertex 
correction. 

Note that, even after the UV renormalization procedure, the Bhabha scattering 
differential cross section still shows poles in {D — 4). This is due to the fact that IR 
divergences are still present and they are regularized by the same parameter D that 
regularizes the UV ones. These divergences are not physical and will be removed 
adding to the process also the corresponding soft real emission. 
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The contributions to the differential cross section appearing in the remainder 
of this work are all UV-renormalized. In the following, whenever discussing the 
contribution of a UV divergent one- or two-loop diagram, we assume as understood 
that the appropriate counterterm graphs have already been added and therefore the 
quantities under study are UV finite. 



4 One-Loop Differential Cross Section 

For completeness, we report in this Section the contribution to the cross section given 
by the interference between the one-loop virtual graphs and the tree-level amplitude. 
The one-loop diagrams contributing to the Bhabha scattering process are of three 
kinds: vacuum polarization diagrams, vertex correction diagrams, and box diagrams. 
Their contribution are discussed separately in the following subsections. 



4.1 One-Loop Vacuum Polarization Corrections 



The two one-loop diagrams containing a vacuum polarization subgraph can be ob- 
tained by inserting a photon self energy correction in the photon propagator, as 
shown in Fig. |21(a) and (b). The photon self energy can be written as 



(23) 



where p indicates the generic 4-momentum entering the self-energy graph. The UV 
renormalized Il^''\p'^) has the following Laurent expansion: 



n, 



(11,0), 



-p 



(24) 



the exphcit expression of the Ill^^'^\—p'^),i = 0, 1 is given in Eqs. ()82I83|) . 

The self-energy diagrams in Fig. |21 (a) and (b) factorize in the product of the 
tree-level amplitude times the function II^q''\—p'^), so that one finds: 



(25) 



The superscript "S"' in the equation above stands for "self-energy" . 

This set of virtual QED corrections is IR finite; therefore, the contribution of 
these diagrams to the differential cross section at order is: 



da^ {s, t, w?) 



a 



1 



st + — + {t- Irn^f 



2Ren, 



(1«,0), 



St + - + (s - 2m^f 



2n, 



+ ^ {{s + tf - 4m^] (Ren(^''°)(s) + J^i'^'\t) ^ 



(26) 



9 



4.2 One-Loop Vertex Corrections 



At the one-loop level, there are four vertex correction diagrams (see Fig. |21(c), (d), 
(e) and (f)). At first, we consider the diagram in Fig|21 (c). The vertex correction 
to the electron current can be written as 



.{11), 



(27) 



where a'^" = —i/2 [y^, 7^^] and t = —{pi —p^Y (space-like) ; the electron and positron 
spinors have been omitted from Eq. (jTTjl . The functions F'^^\t) and F^^\t) are the 
UV renormalized form factors; their Laurent expansion reads 



F2^\t) = Ff''°)(t) + 0(^(Z}-4) 



(28) 
(29) 



The renormalized form factor Fi^''\t) still shows an IR polar divergence in {D — 4); 
the explicit expressions in terms of harmonic polylogarithms are given in Eqs. ()84I85I86|) 
of Appendix 1X1 

The contribution of the diagram in Fig. |21 (c) to the Bhabha differential cross 
section at order can be written as 

daY{s,t) 



a 

(V,{c)) S 



^^/^'^(^,t) + lFf)(M) 



(30) 



where 



2 ' 



i(D-4)(st-2sm2 + s2 



- 2tm^ + 1^) 



ReF[^'^ it) + 2 (^st - Urn" + ^t^^ ReF^ {t) , (31) 



2 St - Asm^ + + -t^ + 4m^ + -{D - 4)t 



+ 2 ( tm'^ + ^t'^ ) ReF^^^\t). 



ReFi^'\t) 



(32) 



Due to the IR pole present in Fi^''\t), the two functions v}^^\s, t),i = 1,2 have the 
Laurent expansions 

V}'\s,t) = + V}''''\s,t) + 0{{D- 4)) . (33) 

The explicit values of the Laurent coefficients can then be expressed in terms of the 
Laurent coefficients of the form factors Fj^^''\t) in Eqs. ()28l29p : 

v}''^-'\s,t) = c,,,{s,t)ReFi'''~'\t) (34) 
V}''''\s,t) = c,i(s,t)ReFf''°)(t)+c,2(5,t)ReFf-^\t) 

+c,,s{s,t)ReFi''''\t) , (35) 
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with 



cii(s,t) = (s + t)2-4m^, (36) 
ci,2(s,t) = ^{st - 2sm^ + - 2tm^ + t"^) , (37) 

ci,3(s,t) = 2(^st-^tm^ + ^t^^ , (38) 

C2,i(s,t) = 2 (^st - Asm^ + + h'^ + im^^ , (39) 

C2,2(s,t) = (40) 
C2,3(s,t) = t(t + 2m2) . (41) 

For convenience of later use in the above equations we write IieFj^^''\t) {i = 
i , / even if the Fl^^\t) are real functions for space-like t. The first term between 
square brackets in Eq. (j3(Jj) comes from the interference between the tree level s- 
channel diagram and the diagram in Fig. |21 (c); the second term originates in the 
interference between the same vertex correction and the tree level t-channel. The 
terms proportional to {D — 4) in Vi(s,t) and V2{s,t) give rise, when multiplied 
by the IR pole in Fi^''\t), to the finite contributions proportional to Q,2(s,t) in 
Eq. (IHHjl . We observe that in Vi{s,t), the term ci,i(s,t) of Eqs. ()34p35j) is identical 
to the coefficient of the s-t-interference term in the tree-level, Eq. (jH}, while the 
term C2,i(s,t) is twice the coefficient of the 1/t^ term in Eq. (jH)). Those relations 
are expected since (t) appears as a multiplicative factor of the tree-level QED 
vertex in the amplitude to be squared. 

The diagram in Fig. |21 (d) can be obtained from the diagram in Fig. |21 (c) with 
the transformations p2 ^ —ps and ^ —pi- Since these transformations leave the 
Mandelstam invariants s and t unchanged, the contribution of the diagram (d) to 
the cross section is equal to the contribution of diagram (c). 

The contributions of diagrams (e) and (f) in Fig. |21 are also identical, since 
the latter can be obtained from the former by exchanging pi ^ — p4 and p2 ^ —Ps, 
which leave s and t unchanged. Moreover, diagram (e) can be obtained from diagram 
(c) by the replacement p2 ^ —p3 which exchanges s ^ t, so that the contribution 
of diagram (e) to the differential cross section at order will be 



V{„ +^ ^,2 

(42) 



daY{s,t) a 



(V,{e)) 



X'(M) + V'^(M) 

3^ St 



According to the definitions in Eqs. (j31l32j) . the expression above involves the func- 
tions F^^''\s) {i = 1, 2); above threshold (s > 4m^) they develop an imaginary part, 
which obviously cancels out in the cross section at this order. In view of that, in 



writing Eqs. ()31l32j) we explicitly indicated that only the real part of the functions 
F^*'^'^ contributes to the cross section. 

The total contribution of the four vertex correction diagrams of Fig. |21 (c)-(f ) to 
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the differential cross section is 
daY{s,t) 



(iiy) 



2"- 

s 



+1 {Vl'\s,t) + V}'^t,s) 



(43) 



4.3 One-Loop Box Corrections 

The one-loop virtual corrections to Bhabha scattering in QED include four box 
diagrams: the direct box diagrams in the t- and s-channel (shown in Fig. |21 (g) and 
(h), respectively), and the corresponding crossed graphs (shown in Fig. |^ (i) and 
(j)). The interference of the diagram in Fig. |21(g) with the tree- level matrix element 
is given by 



daY{s,t) 



dn 



a 



^ReB'^'\s,t) + ^ReBi''\s,t) 



(44) 



where the quantities B^''^ (s, t) and B^''^ {s, t) are functions of the s and t Mandelstam 
invariants (and of the electron mass); they are UV finite, but develop an IR pole in 
{D — 4) (to be canceled out in the physical cross section once the real soft radiation 
is accounted for) so that their Laurent expansion reads 

{s,t)= p^i^f '-'^ is,t) + {s,t) + 0[{D-A)) . (45) 



The explicit expression of the various functions is given by Eqs. ()87I88I89I90|) in 
Appendix El 

The diagram (h) in Fig. Elis the same as diagram (g), in the same figure, once we 
exchange p2 <-> —ps, which is equivalent to the exchange s ^ t. As a consequence, 
the contribution of diagram (h) to the differential cross section at order is: 



dai{s, t) 



dn 



{B,{h)) 



a 
47 



—ReB^i\t,s) + —ReB^^\t,s) 
s t 



(46) 



Interfering the diagram (i) in Fig. El with the tree- level graphs one finds: 
daY{s,t) -2 



dn 



a 
47 



(47) 



where -83^''' {u, t) is another auxiliary function, with the same Laurent expansion 
as Eq. ()45|) . whose explicit expression is given in Eqs. ()91l92j) of Appendix El In 
Eq. (^Tj), the functions B^^\u, t) and B^''\u, t) do not develop an imaginary part in 



the physical region, so that ReB^''\u,t) = B^'^''' {u,t), and Rei?2^'''(M, t) = B!^'^''' {u,t). 

The diagram (j) of Fig.|21can be obtained from diagram (i) of the same figure, by 
exchanging p2 ^ —P3, so that s ^ t; the contribution of this diagram to the cross 



(10, 
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section can then be obtained from Eq. (j47j) by replacing t with s in the arguments 



of the functions B^''^ and B^ 



daY{s,t) 



iB,{j)) 



a 



m 

s 



m 



u, s) ReB^''\u, s) 



(4J 



Summarizing, the contribution of the interference of the 4 diagrams in Fig. |21 (g)- 
(j) with the tree- level matrix element to the Bhabha scattering differential cross 
section is given by: 



daY{s,t) 



dn 



a 
4J 



2 

^ (ReB[''\s, t) + Re5f ^ (t, s) + 5f ^ (m, t) 



>(iO/ 



B^^'\u,t) + ReB\''\u,s) 



(10/ 



(49) 



We checked the expression of the IR pole for the direct t-channel diagram of 
Fig- El (g) against the corresponding result given in |22) where the author calculates 
the one-loop virtual corrections to electron-muon scattering, without neglecting the 
two fermion masses and using a small photon mass A for regularizing the IR diver- 
gences. Setting the mass of the muon equal to the mass of the electron in [2H], we 
find full agreement with our result provided that the dimensional IR pole is replaced 
by 

- (y) • (50) 

We compared also the contribution of the s-channel diagrams to the differential 
cross section (pole and finite part), in the limit m — 0, with the results in [21], 
where the A-regularization is used and the Dirac algebra is in 4 dimensions. If we 
calculate the contribution of the diagrams in Fig. |21 (h) and (j) to the cross section 
using D-continuous dimensions for the Mis but 4 dimensions for the Dirac algebra, 
in the limit m — > we find again agreement with [25, provided the correspondence 
of Eq. (j5(J|) is used. When calculating all the Dirac traces also in D dimensions, 
we obtain, compared to the results in [21], an extra finite term proportional to the 
coefficient of the IR pole in (D — 4). 

Let us recall here that the expression of the one-loop scalar box integral in the 
D-regularization, given in [Hj and used in this paper, is the same as the result of Ap- 
pendix E (b) of [23 (where a small photon mass is still used), once the replacement 
([301) is carried out. 



5 Two-Loop Np = 1 Differential Cross Section 

In this Section we discuss the virtual diagrams contributing to the Bhabha scattering 
differential cross section at order (summed over the spin of the final states and 
averaged over the spin of the initial ones). We limit our analysis to the diagrams 



13 



that include a closed electron loop. Beside the genuine two-loop diagrams, we have 
to consider the interference of the one-loop vacuum polarization graphs with the 
one-loop vertex and box corrections. 



5.1 Two-Loop Vacuum Polarization Corrections 



The two-loop photon self-energy can be written in a way completely analogous to 
the one-loop case: 



(51) 



The explicit expression of the UV-renormalized function nQ^'''(— p^) in terms of 
HPLs is given in Eq. (j93p of Appendix ^ 

The contribution of the two-loop vacuum polarization diagrams in Fig. El (a)-(f ) 
can be obtained by replacing IIq^''' by IIq^''' in Eq. (j^ : 



{2l,S) 



s [ 

1 



st + — + {t- 2m^f 



2Ren;,'')(s) 



+ 



t2 



st + — + [s — 2m 



2\2 



2nf^(t) 



+^ {{s + tf - (Renf )(s) + J^i'\t) 



(52) 



5.2 Two-Loop Vertex Corrections 

The two-loop Np = 1 virtual corrections to Bhabha scattering include the four 
diagrams shown in Fig.|ni(g)-(j) (and the corresponding counter terms), obtained by 
including a vacuum polarization correction in the photon propagator of the one-loop 
vertex diagram. The Lorentz structure of these corrections is identical to the one- 
loop case, Eq. (|77j). with the form factors F'>^^\—p^) {i = 1,2; — = s,t) replaced 
by F^'^^\—p^) {i = 1, 2; —p^ = s, t). The explicit expressions of the UV renormalized 
form factors F-'^''\—p'^) are given in Eqs. ()94I95|) of AppendixEl at variance with the 
one loop case, the two-loop contributions of the graphs in Fig. El (g)-(j) are infrared 
finite, so that no expansion in [D — 4) is needed. 

The full set of two-loop Np = 1 vertex corrections to the squared amplitude can 
be written, similarly to Eq. as 



(i(T^(s, t, rn? 



dn 



2"- 

s 



ivr(M)+^i'f'(^,«) 
+i(vr(.,t)+vf'(M)) 



where the functions V\ and V2'"'' are given by 



(20 



V}^'\s,t) = ci,i(s,t)ReFr'(t) + ci,3(s,t)ReF2^"^(t) 



(20, 



^(20/ 



\/f^(s,t) = C2,i(s,t)ReFr(t) + C2,3(s,t)ReF2^"^(t) 



(20, 



.(20, 



(53) 



(54) 
(55) 
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and the coefficients Ci^i(s,t) etc. are the same as in Eqs. (jnEHH))- 



5.3 Two-Loop Box Corrections 

There are 8 two-loop box diagrams including a closed fermion loop, shown in 
Fig. dSD (k)-(r). 

Let us consider the diagram (1) in Fig. El once we transform the external mo- 
menta according to p4 ^ —pi and p2 ^ —ps, this diagram becomes identical to 
diagram (k). Since this set of transformations leaves unchanged the Mandelstam 
variables s and t, the contributions of diagrams (k) and (1) to the differential cross 
section at order a^lNp = 1) is identical. Similar considerations hold for the pairs 
of diagrams (m)-(n), (o)-(p), and (q)-(r). The diagrams (k), (m), (o), and (q) in 
Fig. El are linked by the same symmetry relations that connect the four one-loop 
box diagrams in Fig. |21(g), (i), (h), and (j), respectively (see the discussion in Sec- 
tion |01)- Consequently, in strict analogy to the one-loop case, the contributions 
of the two-loop box diagrams in Fig. Q to the Bhabha scattering differential cross 
section can be written, similarly to Eq. (|49p . as 



da^is, t, m?) 



dn 



= 

(21,3) As 



2 



(ReSf ^(s, t) + Re^f ^(t, s) + Bf\u, t) 

2 

Re5f V, s)) + ^ (^ReSf ^(s, t) + ReSf '^(t, s) 



Bi^'\u,t) + ReBi^'\u,s)) . (56) 



We choose to factor out in Eq. an overall multiplicative factor —2 accounting 
for the minus sign coming from the closed fermion loops and the identity of the 
contributions from the pairs of diagrams (k)-(l), (m)-(n), (o)-(p), and (q)-(r) of 
Fig. El 

The Laurent expansion of the UV-renormalized functions i?!^'-* reads: 

B?'\s,t) = p^^f '"'H^,^) + Bf^''\s,t) + 0((D - 4)) . (57) 

The expressions of the coefficients of the Laurent expansion, uf^'^"^ (j = —1,0), 
are given m Eqs. (jHEHinU) in Appendix El 



5.4 Two-Loop Reducible Corrections 

The complete set of two-loop diagrams to the Bhabha scattering process containing a 
closed electron loop includes the four 1-loop reducible diagrams, shown in Fig. El (s)- 
(v). They can be obtained from the one-loop vertex diagrams by inserting a one- loop 
photon self-energy correction in the photon propagator connecting the two fermionic 
currents. Since the vacuum polarization corrections amount to a multiplicative 
factor Yi^Q^\—p'^)^ the contribution of such diagrams to the the differential cross 
section at order is easily obtained by appending the appropriate Ii^Q^\—p'^) {p^ = 
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—s, —t) to the various terms in Eq. (j43j) and properly accounting for the {D — 4) 
expansion. In close analogy to Eq. PHj) we write therefore 



{s, t, rn?) 



2 r 



{21, R) 



5^ 



+-{V,\s,t) + V,^{t,s)) 



(5J 



where the V/^ are to be Laurent expanded as in Eq. ()33|) . The coefficients of the 
Laurent expansions are 



(59) 



V!''''\s,t) = c,,{s,t) Re( Ff''°)(t) U^^^''\t) + Fi'^'-'\t) U^^^''\t) ) (60) 

+c,,2(5,t)Re(Ff''-^Ht)n(^''°)(t)) (61) 
+Q,3(s,t)Re(Ff'°)(t)n(^''°)(t)) , (62) 



where the coefficients Ci_j{s,t) are given by Eqs. 
nj)^'^(t) and of the F^^^^ {t) by the Eqs. 



the Laurent coefficients of 



5.5 Products of Two One-Loop Corrections 

Finally, to conclude our catalog of contributions to the Bhabha scattering cross sec- 
tion at order a^lNp = 1), we need to discuss the interference between the one-loop 
vacuum polarization diagrams of Fig. |21 (a), (b) and the one- loop vertex and box 
graphs of Fig. |21 (c)-(f) and (g)-(j). Since the one-loop vacuum polarization in- 
sertion factorizes in the tree-level diagram and the multiplicative factor Il^''\—p'^) 
(— = s,t), the contribution of these interferences among one-loop diagrams is eas- 
ily obtained by replacing the tree-level photon propagators 1/s and 1/t in Eqs. f)lH| 
and (gH) by Ren[,^'^(s)/s and Il^Q^\t)/t, respectively. 

In this way one obtains for the interference of the self mass graphs in Fig. (j2I) (a), 
(b) and the vertex graphs in Fig. Q (c)-(f) 



(i(T^(s, t, m?) 



dn 



(S,V) 



a 

2— Re 

s 



lyf)(M)(nr(. 
+^^fH^,t)(nr(t) 

+lv/^')(M)(nr(t) 



(63) 



where the functions v}^''\s,t) are the same as in Eqs. 



. As they have the 



Laurent expansion of Eq. the above cross section can be Laurent expanded in 
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the same way. One has for instance, in the notation of Eqs. (j33|) and (j24j) : 

+vr\t,s)(ur\s)y+v<'^'-'\t,s)(ur\^ 

+0((D-4)), (64) 



and similar expressions for all the other terms of Eq. ()63|) . 

For the interference of the self mass graphs in Fig. El (a), (b) and the box graphs 
in Fig. 121 (g)-(j) one finds, in analogy with Eq. 



(icr^(s, t, m^) 



(S,B) 



a . 
47" 



s V 

+5f(.,t)-B?')K.))(nrO 

-Bi-\u,t) + Bi-\u,s)){ni-\t)y 



(65) 



where the functions -Bj-^'^(s, t) are the same as in Eqs. ()44I47|) . with the Laurent 
expansion of Eq. A similar expansion holds for Eq. ()65|) : in analogy with 

Eq. ()64|) . one has for instance: 



B['\s,t){n^^\s)y = ^^B[--^\s,t){ui 

+Br\s,t)[nr\s)y+B['^'-'\s,t){nr\s) 

+0((D-4)), (66) 
and similar expressions for all the other terms of Eq. (jHSj)- 



6 Conclusions 

In this paper, we considered the virtual radiative corrections to the Bhabha scatter- 
ing differential cross section (summed over the spins of the final states and averaged 
over the spins of the initial ones) in pure QED; we calculated the contribution of all 
the one-loop diagrams (order a^) and of the subset of two-loop diagrams involving 
an electron loop (order a'^{Np = 1)). 

The provided result does not rely on any kind of approximation; the full depen- 
dence of the cross section on the electron mass m and on the Mandelstam invariants 
s and t was retained. 

Both UV and IR divergences were regularized within the framework of dimen- 
sional regularization. The UV renormalization was performed in the on-shell renor- 
malization scheme. 
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The reduction of the integrals to the Master Integrals was carried out employing 
the Laporta algorithm. The Mis involved in the calculation of the cross section 
were evaluated in two previous works [HI HH] , by means of the differential equation 
method. 

All the results are written in terms of 1- and 2-dimensional HPLs of maximum 
weight 3. 

The differential cross sections at the one- and two-loop level can be obtained 
adding up the various contributions discussed in the paper. 

In particular, the one-loop virtual cross section is the sum of the contributions 
given in Eqs. ()26f43p49|l respectively: 



dn 



da^is, t, m? 



dn 



(iiy) 



dermis, t, m? 



dn 



; (67) 



the two-loop differential cross section is the sum of the expressions given in Eqs. 

EHEiinnEsi): 



(i(T^(s, t, rn?) 

dn 



+ 



dn 



{21, s) dn 

(icr^(s, t, m?) 



dn 



+ 



dn 



dn 

d(T2{s, t, m? 



{2l,B) 



iS,V) 



dn 



(S,B) 



(6J 



The differential cross section still includes single poles in {D — 4), of IR origin, 
to be canceled by the inclusion of the contribution of the diagrams involving the 
emission of a real soft photon, not discussed in this paper. 
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A Auxiliary Functions 



In this Appendix, we give the exphcit expressions of the auxihary functions used 
throughout the paper. They involve 1- and 2- dimensional Harmonic PolyLogarithms 
(HPLs), introduced and discussed in [1^1 CHI HZI CHI HH E] , of arguments x, y related 
to the kinematical Mandelstam variables s, t by the relations 



-t = (p^ - = Q = m , y = " >^j= . (70) 



y ' v^g2 + 4^2 + ^ 

For completeness, we introduce also 

. ^2 ,.2 2 (1-^)' + 4^2 _ 7^2 

- M = (pi - P4) = V" = m ; , z = — == ; (71) 

the variable 2; is indeed related to x and y through the relation 

- M = s + t - W = m -(x + y) x + - = m^- — . (72) 

X \ yj z 

If s = —P^ is space-like (s negative, positive), x as given by Eq. (jU^ is real and 
positive, varying from x = lats = Otox = Oats = — 00, and similarly for the 
variable y. 

In the explicit formulas which will follow, the one-dimensional HPLs are written, 
according to 15 , as H{0, 1; x), H{—1, 0; y) etc.; examples of two-dimensional HPLs 
are G{—l/y;x),G{—y,0,0;x) etc., for their definition see [6^. 

The arguments of the auxiliary functions will be taken to be space-like, (— P^ in 
the case of a single argument, (— P^, —Q^) if the functions depend on two variables) 
so that the corresponding x and y vary both in the (0, 1) range and all HPLs are 
real. The time-like regions can be recovered by analytic continuation. 

Let us consider, for definiteness, the variable s = —P^. The continuation to the 
time-like region is performed, according to the usual ie-prescription, by setting 

p^ = -s-ie, (73) 

when s > 0. For < s < 4m2, a/^^^s"^-^ = —i^/s and Eq. gives x = r, where r 
is the phase factor 



V4m2 - S + i^/s _ ^^20 



V4m2 — s — i\/s 



e'"f' , (74) 



with 



arctan 4 / . (75) 

\/ 4^2 _ g 



Above threshold, s > 4m^, we define 



x' = , (76) 

\/s + \/ s — Am? 
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with x' = 1 at s = Aim? and x' = at s = cxd and the continuation in x is performed 
by the replacement 

X = —x' + ie . (77) 

Remembering that 

if(0; x) =liax , 

the continuation from space-hke s to s above threshold is then performed by 

H{0; x) = H{0; -x' + ie) = H{0; x') + in . (78) 

Another less trivial example is 

H{-l,0;x) = H{-l;x)H{0;x) - H{0,-l;x) 

= H{-1; -x' + ie)H{0; -x' + te) - H{0, -1; -x' + ie) 

= -H{l;x')(^H{0;x') +iiT^ + H{0,l;x') , (79) 

where an algebraic property of the 1-dimensional HPLs (shuffle algebra) has been 
used in the first line and the relation reversing the sign of the argument in the last 
line. For more details on the analytic continuation of the HPLs we refer to the 
already quoted papers. 

The contribution of the two-loop box diagrams to the differential cross section 
includes a subset of the 2-dimensional HPLs introduced in |6|; this subset can be 
expressed in terms of Nielsen polylogarithms (whose argument depends however on 
both X and y, i.e. on both s and t), exphcitly given in Eqs. (C.21,C.22,C.27-C.30) 
of jB] . As an example of analytic continuation which applies to those cases, consider 
for instance Li2{—x/y); if x,y correspond to s,t when they are both spacelike, the 
continuation to physical (timelike) s and spacelike t implies x —x' + ie, with 
x' < y, so that one has simply Li2{—x/y) Li2{x'/y). The case of timelike t, 
spacelike s, on the other hand, implies y — y' + ie, with y' < x, so that the 
continuation is given by 

Similarly for Lis we have: 

x\ T- ■\ T- fy'\ T 3 fy'\ oA/r,\ 1 f y' 



U,{ — ]^U,(- + ie\ = LigM^ +-ln^M^ -2C(2)ln 

\ yj \y J \xj 6 \xj \x 

In the following subsections we provide the explicit expressions of the auxiliary 
functions. 
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A.l Auxiliary Functions for the One-Loop Cross Section 



n, 



(l/,0). 



4 4 
+ 



9 3(1 -x)2 3(1 -a;) 



1 4 

+ 



3 3(1 



n, 



-p2) 



(1 

14 



16 



16 



27 9(1 9(1 -x) 



1 1 
f- 



+ 



16 



+ 



3(1 (1 -x)2 6 
1 



(82) 
C(2) 



18 9{l-xf 3(l-a;)2 3(1 - x) 
12 1 " 



H{0;x) 



+ 



6 3(1-,t)3 (l-a;)2_ 
14 2 



3 3(l-a;)3 {1 - x)\ 



H{0,0;x) 
Hi-l,0;x)} , 



(83) 



(11,0). 



,{ll,0) 



1 - 



1 - 



1 



1 



-1-1 

2 



1-- 



(1 - x) (1 + x) 
1 1 1 

H{0;x) 



HiO;x) , 



2 {l+x)i 
1 



1 



(84) 



2L (l-x) (l+a;)J 
-2if(-l,0;a;)], 
1 

"2 



(l-x) (l + x) 



[C(2)-2i/(0;x)-i/(0,0;a;) 



i/(0;x) . 



(85) 
(86) 



-48 



8 



32 



32 



x'^{l — y)^ x^{l—y) x{l — y)^ x{l — y) 



16 32,x 32a; 8x 
+ - + 16x4 



5X 



X {1-yy {1-y) 



(l-l/)2 (l-y) 



8 



8y 



?/(l + x) ?/(l-x) y (l + x) (l-x) 

96 96 80 

■(l + a;)(l-t/)2 + (l + a;)(l-y)^(l + x) 



+ 8?/ 



32 



32 



16 



(l-x)(l-y)2 (l-a;)(l-|/) (l-x) 



32 



32 



{1-yy (l-y) J 



H{0;x), 



(87) 



S^VP^-Q^) = C(2) 



32 



x'^{l + y) x^{l-y) x'(l + y)-^ x(l + y)^ 
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88 40 12 32a; 48a; 

x{l + y) x{l-y) ^ {l + yy (i + y)^ 

(1 + y) (1-y) (1 + y) (1 - y) 

4 6 ^iQy^ 32 

y{l-x) y {l-x) (1 + 2^) 

64 96 208 48 \ 

+ (rTW" (TT^+lTT^" (1^" ) 

/ 4 4 8 8 2 

\ x'^{l — yY x'^{l — y) x{l — yY x{l — y) x 

8x 8a; „ 4a;^ 4a;^ 

1 I- 2r -I 

(l_^)2 + (i_^)+ ^ (l-yY {l-y) 

4 _ 8 ^2^4y_8y 

y{l + x) y{l-x) y (1 + x) (l-x) 

32 32 16 

+ ^ ~ (1 + X)(l - + (1 + x)(l - t/) + (T^ 

~ (T^) ^) + ('^ - ^(TT^ + ^(iT^ 
"x" (TT^+ (TT^" ^^(1+^ 

4 16x 16a; a ^ a 

(rT^ + (rT^" ""'y' ^+(1+^ 

(i + y)y V ^ y 

8u 64 \ , , / 4 



4 



8y 64 \ / 4 

4 16 16 4 16x 

^(1-1/) x{l-yY x{l-y) X [l-yf 

16a; , 4a;2 4a;2 4 2 

/^x 1 1 

(1-y) {l-yf (1-y) y(l + a;) y 

4y 48 48 

■(1 + a;) ^+(l + a;)(l-y)2 (l + a;)(l-y) 

8 16 16 « 



"(1 



(1 + a;) (1 -x)(l -?/)2 ^ (1 -x)(l -I/) (1-a;) 

(r^ + (T^)«(»;-W'';«) + (i«+?(i^ 

8 32 32 8 32a; 

'(1-y) a;(l-y)2 a;(l - y) a; (1 - y)^ 



X- 



22 



32x 8a;2 Sx^ 8 
- 8a; - -r + + 



(l-y) {\-yf (1-2/) y{\^x) y 

8y ^ 96 96 



(1 + x) " (l + x){l-yy {l + x)(l-y) 

16 32 32 16 

+ 



(1 + x) (l-x)(l-y)2 (l-x){l-y) (l-x) 

32 32 \ / 4 

+ - -)H{0;x)H{l;y)+ 16 + - - 4x 



(l-y)2 (1-2/); ^ ' ^ V ■ ^ 
+ RI^-^(I^-^^-(If^)^^°'°^^^ 



x2(l + 2/) a;2(l-2/) a;(l + 1/)^ a;(l + 2/)2 

22 10 4 8x 12x 

+ T H t:^ ^ + 



x(l + y) x(l-2/) X (l + y)3 (l + y)2 

22x lOx , 2^2 2x2 2 ^ 
+ T. ^ + 4:X + - h 2y 



(1 + y) (1-y) (1 + y) (1 - y) y 

16 24 52 12 \ 

+ 1-32+ ^ 4 16 24 



x'^{l + y) x'^{l—y) x{l + y)'^ x{l + yy 
44 20 8 16x 24x 

+ -7Z T H -7Z + 



x{l + y) x{l-y) X (l + 2/)3 {1 + yf 

Ux 20x ^ 4^2 4^2 4 ^ 
+ 7 T + Sx + -. ^ - -. ^ h 4^/ 



(1 + y) (1-y) (1 + y) (1 - y) y 

32 48 104 _ 24 \ 



5(i'.-i)(_P^,_g^) = ( -48- -- + 16x4 



16 16 16 _ 16x- 



2 



x2(l — y)2 x^(l—y) X (1— y)2 

16x2 8 8 8 8y 



Z/(l + a:^) z/(i-2^) y (l + x) 
8y 64 64 

+ 8y - 7r—T7^-T, + 



(l-x) (l + x)(l-2/)2 (l + a;)(l-2/) 

80 64 64 16 

'(1+x) " (l-x)(l-l/)2 + (l-x)(l-|/) + (r^ 

64 64 \ 



fuo). 9 ^/ 16 16 64 96 



x2(l + 2/) x2(l — I/) x(l + 1/)^ x{l + yy 
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80 16 12 64a; 96a; 
+ 



x{l + y) x{l-y) X {l + yy (1 + 1/) 

SOx 16x 16^2 16x2 

+ ^ + 20a; + 



2 
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2 6 2y 2^/ 20 

lOy + 



y(l-a;) y {1 + x) (1 - x) ^ + x) 
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-I 1 1 qo 
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\ X y{l — x) y (1 — a;) 

16 \^^,^ , 16 16 2 

+7i N 2;) + 12 - ,^ , + 



(1 — a;) y ' \ a;(l + y^ x{l + y) x 

16x 16x ^2 32 
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32 \ , f^, 32 32 4 

.——-\H{0;y +[24- + -- 
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24 



, 4 ^ 2 2 2 2y 

+ ( 12 + --4X+— r + ^- r - - + , 

X 1/(1 + a;) 1/(1 — 1/ [1 + X) 



2?/ 



(1-x) 
+ ( -20 + 



-2|/- 



20 



(l + x) (l-x) 
4 4 



H{0,0;x) 

16 24 

+ 



x2(l + y) x^(l-y) x(l + y)^ x{l + y)^ 

20 4 4 16a; 24a; 20a; 
r + 7^ + 



x{l + y) x{l-y) X (l + y)3 (l + y)2 (l + y) 

Ax^ 2 32 



4a; 4a;^ 
+ TZ T + 4a; + 



48 



(l + y) (l-y) y 
64 8 \ 



2 

- + 2|/ 



(l + y)2 



(l + l/)2 (l + y) (1-y); 
8 



+ -40 + 



H{OAy) 

32 



+ 



48 



a;2(l + |/) a;2(l-|/) x(l + |/)3 a;(l + |/)2 

40 8 8 32a; 48a; 40a; 

^— r + — — - + 



x{l + y) x{l-y) X {1 + yy (l + y)' (l + y) 



+ 



ix 



ix 



ix 



(1-y) 

96 



5a; 



(l+y) (1-y) y 
128 16 \ 



4 ^ 
- + 4y 



64 



(l + y)2 (l + y) (1-y); 



iJ(0,l;y), (90) 



32 



32 



^2(1- y)2 a;2(l— y) a;(l — y)^ a;(l — y) 



32a; 



+ 



32a; 



8a;2 



+ 



8a;2 



(l-y)2 (1-y) (l-y)2 (1-y) 

32 32 

+ rr. r4 



96 



(l + a;)(l-y)2 (l + a;)(l-y) (l-a;)(l-y)2 
96 32 32 



(l-x)(l-y) (l-y)2 (1-y) 



H(Q-x), (91) 



= C(2) 



+ 



+ 



32 



48 



x2(l + y) x2(l — y) x(l + y)3 a;(l + y)^ 

48a; 8a; 



24 



x(l + y) 


a;(l- 


24a; 


8a;2 


' (1-y) 


(1 + 1/) 


2y 


2y 


' (l + a;) 


(1-a;) 


96 


48 



+ 



32x 



5a;2 2 2 

+ 7^ 7 + 



12 



64 



16 



(l + y)2 (l + y) (1-y) 



12 



+ 



4 



a;2(l — y)2 



25 



— y) x(l — yy x{l — y) x (1 — y)' 



+V-. r - 4x - + r + 



8y 



(l + x) 



16 



(1-y) + y 

32 32 



(l+x) 



f7T^li^(0;x)+(-8 + 



+ 



16 



+ 



x{l+yy x{l + y) 
16 \ 

my) 



(l+2/)2 (1+1/) (l+y)2 (1+y) 

16 16 16x 16a; 

+ 



+ -16 + 



x{l + yy x{l+y) {1 + yy {1 + y) 



32 32 \ , / 4 



{1 + yy (l + y) 



4 4y Ay 24 



y(l-x) (l+x) (1-x) (l+x) (1-x) 



y{l + x) 

H{-l,0;x) 



4 4 

+ 



16 16 

+ 



+ 



x^{l — y)"^ x'^{l — y) x{l — y)"^ x{l — y) 
16x 16x 4x2 2 



(l-y)2 (1-y) (l-|/)2 (1-y) y(l + x) 
2 2|/ 2y 16 



y(l-x) (l + x) (1-x) (l + x)(l-|/)2 
16 12 48 48 



+ 



+ 



(l+x)(l-y) (l+x) {l-x){l-yy (l-x)(l-y) 

li/(0;x)if(0;i/)+ ( -8 



4 16 

+ 



(1-x) {1-yy (1-1/) 



32 



32 



+ 



32x 



x2(l— |/)2 x2(l — y) x(l — |/)2 x(l — y) (1 — y)2 



32x 



8x2 



<iX 



(l-y) (l-2/)2 (l-y) 1/(1 + x) 2/(1 -x) 



4y 



4y 



+ 



32 



32 



(l + x) (1-x) (l + x)(l-?/)2 (l + x)(l-y) 
24 96 96 8 

^TZ T - TZ rjZ TTT + 



(l + x) (l-x)(l-y)2 {l-x){l-y) (1-x) 



32 



(l-y)' (i-y) 



32 \ / 2 

'H{0;x)H{l;y)+i4- 



2y 2?/ 



?/(l + x) 

12 4 

fT: rlii'(0,0;x) 



7/(1 -X) (l+x) (1-x) (l+x) (1-x) 



+ 4 



+ 



12 



x2(l + 2/) x^{l-y) x(l + 2/)3 x(l + 2/)2 



26 



6 



8x 



12x 



2x 



x{l + y) x{l-y) {l + yf {l + y^ {1 + y) 



+ 



6x 



2x^ 

+ 



16 24 

+ 



(1-y) (1 + y) (1-y) + (l + y)^ 

12 



(1 + y) {l-y)J 



^ ]H{0,Q;y)+(8 ^ 



+ 



16 



24 



x^{l + y) 

4 12 

+ 



+ 



x^(l — y) x{l + yy x{l + yy x{l + y) x{l — y) 

16x 24x Ax 12x 4,x^ 

+ 



+ 



(l+|/)3 (l + |/)2 (1+y) (l-y) (1 + y) 

4a;2 32 48 24 

+ 



(1-y) (l + y)3 (l + y)2 (1 + y) 
i^(0,l;y) . 



(l-y) 



(92) 



A. 2 Auxiliary Functions for the Two-Loop Cross Section 



13 



13 



24 6(l-a;)2 6(1 - x) 



+ 



(l-a;)4 {1-x 



(l-x) 
^ 4 



4 (l-a;)3 2(l-a;)2 (l-x) 



i/(0;j;) 



16 



+ 



3 3(1 -x)3 (l-a;)2 
26 23 2 



ii'(-l,0;x)+ 2 



3(1 -a;)4 



3(1 -x)3 3(1 -x)2 3(1 -x) 
^ V(l.O; • " 



,2 8 

+ - + 



(1 - X)2 J 

16 

■3(1 - ,X)2 

16 



+ 



3 3(1 -x)3 
32 



3 3(1 -x)4 3(1 -x)3 

2 



i/(0,-l,0;x) -//(-l,0,0;a;) 



3(l-a;)3 3(1 -a;) 



3 3(1 -x)4 
i/(l, 0, 0; x) - H{0, 1, 0; x)] , (93) 



383 



+ 



49 



108 ■ 9(l+a;)2 9(1 + a; 
229 11 



49 /I 9^ 



196 



x] 



6(l + a;)2 


2(1 + a;)/ 


563 


59 


' 108(1 + x) 


27(1 -x) 


163 

+ -— 





C(2) + 



265 



3(l + a;)3 
89 



216 9(1 



19 



- X ) 
62 



6(l + a;)2 
124 



18(1 + x)2 6(1 + a;) 



^(0;x)+(gg + g^^^^^^ 9(l + a;)3 

H(0,0;x)+{ -T + 7 -r 

' ^ ' ' ^ * 6 (l + x)5 (l + x)4 



27 



11 



+ 



+ 



+ 



(l + a;)3 2(l + a;)2 12(1 + a;) 12(1 - a;) 
+H{Q,Q,Q;x)\ , 

17 / 34 68 



C(2)//(0;x) 



(94) 



17 



+ - 



+ 



33 



3(l + j;)2 3(1 + x) V C^ + x)^ (l + x)3 (l + x)2 
1 \ / 31 31 59 



(1 + x) 
17 



C(2) + 
H{0;x) + 



+ 



3(l + a;)3 2(l + a;)2 18(1 + a;) 

22 44 23 



+ 



9(l-a;)y V > y ■ ^ 3(l + a;)4 3(l + x)3 3(l + a;)2 

6 15 21 



3(1 



^)i.(0,0;.)+( 



+ 



4(l + a;)2 8(1 + a;) 8(1 - x) 



{1 + xf (l + a:)4 2(l + a:)3 
^ ^C(2)i/(0;a;) + i^(0,0,0;a:) , (95) 



16 



+ 



32 



28 



56 _ 

y ~ 3x2(1-^)4 ' 3x2(1-^)3 9x2(1-^)2 

20 64 128 16 



9a;2(l-y) 


' 3x(l 


-y)4 3a;(l-y)3 ' 


9x(l-y)2 


176 


40 


64a; 128a; 


16x 


' 9a;(l-|/) 


9a; 


3(1_^)4 1 3(l_^)3 


9(l-|/)2 


176a; 


40a; 


16a;2 32a;2 


28a;2 


9(l-y) 


9 + 


3{l-yy 3{l-yf 


' 9(l-y)2 


20^2 

+ r + 


20 


20 20 


20y 



9(l-y) 9y(l-a;) 9y(a; + 1) 9y 9(1 - x) 
20y 20y 64 



128 



9(a; + l) 9 3(1 - a;)(l-y)4 3(1 - a;)(l-y)3 

112 80 88 

^— — - + 



9(l-a;)(l-y)2 9(1 -a;)(l-2/) 9(1 - a;) 

64 64 128 

+ 



+ 



(l-t/)4(a; + l) 3(l-y)4 (l-y)3(a; + i) 
128 16 80 



3(l-y)3 3(l-y)2(x + l) 9(l-y)2 
208 272 248 



+ 



3(l-y)(a; + l) 9(l-y) 9(x + 1) 

16 40 8 

4 



H{0;x) 



3a;2(l-y)5 3x2(l-y)4 x^{l-y)^ 3x^{l-y)^ 

4 64 160 64 
+7r^. 7 + - T^^z - + 



3a;2(l-y) 3x(l-y)5 3x{l~yY 3x{l-yf 
^ 64 16 8 + 



3a;(l-y)2 3a;(l-y) 3a; 3(l-y)5 3(l-y)4 



28 



64 



64 



16 



8x 



3(l-y)3 3(l-|/)2 3(l-y) 3 3(l-2/)5 



40x2 



4x2 



4x2 



3(l-y)4 3(l-y)2 3(1-1/) 3i/(l-x) 



1 



4 



4y 



4y , 4y 



3y(x+l) 3y 3(1 - x) 3(x + 1) 3 

64 160 32 
rr. - 777^ rr. rT + 



3(l-x)(l-|/)5 3(l-x)(l-y)4 (l-x){l-yf 
16 32 8 

+ 777^ 77^ TTT " 777^ 77^ 74 



3(l-x)(l-2/)2 3(l-x)(l-2/) 3(1 -x) 

64 64 160 160 

f:77:^ 7t4 



(l-y)5(x + l) 3(l-z/)5 (l-y)4(^ + l) 3(l-y)4 
160 32 80 112 

f7777 7^ - 77 7777 7T + 



3(l-y)3(x + l) 3(l-y)3 {l-yy(x + l) 3{l-yy 
32 40 



3(l-y)(x+l) 3(x+l) 



HiO;x)H(0;y) 



(96) 



9 45x2(1 -y)4 45x2(1 -y)3 135x2(1 -y) 

152 64 128 512 

+ 



2 



135x2(1-7/) 45x(l-|/)4 45x(l-y)3 135x(l-y) 
704 304 64x 128x 

+ T77^^7 7 + T77^ + 



2 



135x(l-|/) 135x 4:5(1 -y)^ 4:5(1 - y) 
512x 704x 304x 16x2 

+ T77377 7 + 



135(1 -y)2 ' 135(1 -y) ' 135 45(1 - y)^ 
32x2 104x2 152x2 56 

+ T377 7^ + 



45(1 -y)3 135(1 -i/)2 135(1 -y) 9y(x + 1) 
56 152 56y 56y 152y 

^9y(x + l) ^ 13% ~ 9(x + 1)2 ^ 9(x + 1) ^ ~135 

32 64 32 

+ 



2 



45(1 -7/)4 45(1 -|/)3 3(1 -|/)2(x + 1)2 
32 224 32 

+ 7777^ .n. ... +77^ 777 + 



2 



3(l-z/)2(x + l) 27(1 -y)2 3(l-y)(x+l) 

32 1216 16 16 

+ 



3(l-y)(x + l) 135(1 -y) 3(x + l)2 3(x + l) 
,'64 32 64 56 

+ ^-7r^777 7T + 



3 3x2(1 -7/)4 3x2(1 -|/)3 9x2(1 -|/)2 

40 176 352 140 



9x2(1-1/) 3x(l-y)4 3x(l-y)3 9x(l-y)2 
388 16 16 14 112x 

+ 7r^. 7 + „ . . - „ . . .X + 7^ + 



9x(l-y) 3x(y + l)2 3x(y + 1) 9x 3(1 - y)^ 



29 



224a; 124a; 212a; 16a; 16a; 

3(1 - yf + 9(1 - yf + 9(1^ + ?>{y + 1)^ %TT) 

26a; 16a;2 32a;2 28a;2 20a;2 

+ 



9 3(1 3(1 9(1 - y)^ 9(1 - y) 
10 112 224 116 

^9y(l-a;) ~ 3?/(x + 1)^ ^ 3y(a; + l)3 ~ 3y(x + l) 
22 2 lOy 112y , 224y 



2 



9y(a; + l) 9y 9(1 - a;) 3(a; + 1)^ 3(a; + 1)^ 
llQy ^ 22y ^ 2|/ ^ 32 



4 



3(a; + l)2 9(a;+l) 9 3(l-x)(l-|/) 

64 56 40 
+ rvz + 



3(1 - a;)(l - yf 9(1 - a;)(l - y^ 9(1 - a;)(l - y) 

44 192 160 304 

+ 



9(1 -a;) (1 - t/)4(a; + 1)2 (1 - y)4(a; + i) 3(1 - y)^ 
384 320 608 

+ 



3 



(l-|/)3(a; + l)2 (l-^)3(a; + l) 3(1 - 2/)= 

64 128 368 

+ 



(1 - y)2(a; + 1)4 {l-y)^{x + l)^ 3(1 - y)2(a; + 1) 
184 244 64 

+ .n. . .. -777-. 77; + 



2 



3(l-y)2(a; + l) 9(1 - y^ (l-y)(x+l)^ 

128 208 296 

+ 



(l-y)(a;+l)3 3(1 -y)(a;+ 1)2 3{l-y){x+l) 

668 32 64 344 868 

+ . .... . . . 



9(1 -y) (a; + 1)4 (x + l)3 3(x + l)2 9(a; + 1) 
32 32 / 124 40 



3(y + l)2 3(y + l);"^^ V 3 3a;2(l-y)4 

80 278 82 176 

+ 77^r^. ^ + 



3a;2(l-y)3 27a;2(l-|/)2 27x2(1 ~ y) 3a;(l - y)^ 
352 452 1132 242 

+ 



3x{l-y)-3 27,T(l-y)2 27a;(l - y) 27x 

176a; 352x 452a; 1132a; 242a; 

' 3(1 - yy " 3(1 - y f + 27(1 - y)^ + 27(1 - y) + ^ 

40a:2 80a;2 278a;2 82a;2 

+ 



3(1 -y)4 3(1 -y)3 27(1 - y)2 27(1 - y) 

70 112 56 358 

+ 



27y(l-x) 9y(a;+l)3 3y(x + 1)2 27y(a; + 1) 
^130 70y 112y ^ 56y 358y 



27y 27(1 -a;) 9 (a; + 1)3 3(a; + 1)2 27(a; + 1) 
130y 64 128 



27 (l-x)(l-y)4 (l-a;)(l-y)3 

1664 64 464 

+ 



27(1 -a;)(l -y)2 27(1 - a;)(l - y) 27(1 - x) 



30 



512 



160 



1024 



Sil-y^x + l) 3{l-yy 3(1 - + 1) 
320 64 32 



3{l-yf 3{l-yy{x+lf {l-yy{x + iy 
16 616 64 



3(1 + 1) 27(1 3(1 -I/) (a; + 1)3 
32 496 2056 



{l-y){x+iy 3(l-y)(x+l) 27(1-1/) 

32 16_ 1912 \ /46 

■3(x + 1)3 {x + 1)2 + 27(a; + I) J ^ ' '^'^ ^ 



16 



+ 



8 



8 



+ 



45x2(1 -y)5 9x2(1 -y)4 9x2(1 -y)2 9x2(1 -y) 

2 Ay 64 32 32 



+ + 



15x2 15x 45x(l-7/)5 9x(l-t/)4 45x(l - 1/)3 
208 16 8 148 Axy 



+ 



45x(l - 1/)2 15x(l - y) 3x{y + 1) 45x 15 

64x 32x 32x 208x 

+ 



45(1 -?/)5 9(l-?/)4 45(1 -?/)3 45(1 -y)2 

148x 



16x 8x 

+ T^. T + 



16x2 8x2 

+ 



15(1 -y) 3(y + l) 45 45(1 - y)^ 9(1 - y)^ 



8x^ 2x^ 2x2 2 46y 2y^ 



9(1 -y)2 9(1-1/) 15 y 15 15 45(1 -y)^ 
16 128 656 52 



9(1 -y)4 15(1 -y)3 45(1 - y)2 5(1 -y) 



16 



3(Z/ + 1) 
4 



H{0;y)+(-^4- 



+ 



20 



3 3x2(1 -y)5 3x2(1 -y)4 

2 2 32 
+ ^r-^. r + 



x2(l-y)3 3x2(1 -y)2 3x2(1 -y) 3x(l - y)^ 

80 32 32 8 8 

3x(l — yY 3x(l — y)3 3x(l — y)2 3x(l — y) 3x 
32x 80x 32x 32x 8x 



+ 



+ 



3(1 -y)5 3(1 -y)4 3(1 -y)3 3(1 - y)2 3(1 - y) 

8x 8x2 20x2 
+^ + 



4x2 2x2 

+ 



3 3(1 -y)5 3(1 -y)4 (1 - y)3 3(1 - y)2 
2x2 2^2 2y 2y 



3(1 -y) 3y(l-x) 3y(x + 1) (1 - x) 3(x + 1) 

32 80 16 

t 



4y 



3 3(l-x)(l-y)5 3(l-x)(l-y)4 (l-x)(l-y)3 
8 16 4 



3(l-x)(l-y)2 3(i-x){l-y) 3(1 - x) 

32 32 80 80 

^— + 



(l-y)5(x + l) 3(1 -y)5 (l-y)4(a; + i) 3(1 - y)^ 



31 



80 16 40 56 

3(1 - yf{x + 1) +3(1 - yf " (1 - yfix + 1) + 3(1 - yf 
16 52 ,,,, , / 16 

32 80 16 8 

+ 



3x2(1 3x2(1 x2(l-|/)3 3x2(1 -y)2 
4 4 128 320 

+ 



x2(l — y) 3x2(7/ + 1) 3x(l — yY 3x(l — y) 
128 128 52 16 



4 



3x(l — |/)3 3x(l — |/)2 3x(l — I/) 3x(|/ + 1) 
8 44 8 128x 320x 

+ -77-, TV 7771 + 



3 



x(y + l)2 3x(y + l) X ?,{l-yf 3(1 -y) 

128x 128x 4x 16,T 

+ 



4 



3(1 3(1 -y)2 {1-y) 3(1/ +1) 



3 



^x 44x 8x 32x2 80x2 

+ TTT -T + ^ + 



(?/ + l)2 3(7/ + 1) 3 2>{l-yf 3(1-1/) 
16x2 8x^ 4x^ 4x^ 

+ 7: + 7^. 777 - 7777 7 " 77^^ + 



4 



(l-2/)3 3(l-2/)2 3(1-2/) 3(2/ +1) 32/(1 -x) 
8 4 8y 8y 4y 

+ 3y(x + l) ~ % ~ 3(1 -x) ~ 3(x + 1) + "s" 

128 320 64 

+ 7777 777 7^- 7777 777 77 + 



3(1 - x)(l - yf 3(1 - x)(l - yY (1 - x)(l - yf 
32 64 16 

^3(l-x)(l-2/)2 ~ 3(l-x)(l-2/) + 3(1 -x) 
128 128 320 

+ 7777 77 + 



(l_,^)5(a, + l) 3(l-?/)5 (l-y)4(,x+l) 

320 320 64 160 

3(1-2/)^ ~ 3(l-2/)3(x+l) + 3(1-2/)^ ~ (1-^)2(^+1) 
224 64 8 80 

+ 



3(l-2/)2 3(l-2/)(x + l) (1-2/) 3(x+l) 
32 16 104 \ , /32 2 



3(2/ + 1)3 (2/ + 1)2 3(2/ + 1)7 ' ' ' ' V15 15x2 
4 4v 16 4x 4xv 16x 2x2 2 

+ T^ + T^- — + — + ^- — + — + 



15x2/ 15x 15x 152/ 15 15 15 Iby 
16 16y 2y2\ /112 56 



2 



^5 



15y 15 15 y ' V 9 3x2(1 -y) 

140 32 4 8 

+ 



3x2(1 -y)4 a;2(l-2/)^ 3x2(1 -y)2 3x2(1 - y) 
80 200 896 152 



x(l-y)5 x{l-yf 9x{l-yf ?>x{l - yf 
68 34 80x 200x 896x 

+ 7^^7 7 + 7^ + 



3x(l-y) 9x {l-yf {I - yf 9(1-2/)= 



32 



152,T 68a: 34a; 56x^ UOx^ 
3(l-i/)2 3(1-1/) ~ ~9~ " 3(1 + 3(l-7/)4 

32a;2 4a;2 Sx^ 4 4 

+ 777^ + 



(l-y)3 3(l-?/)2 3(1 -y) 9y(l-,T) i/(a;+l) 
^20^ ly , ly 20y 25(3 



9y 9(1 -x) (x + l) 9 3(l-a;)(l-y) 
640 1408 64 

+ 777^ 77^ 77 - 777^ 77^ 7^ + 



5 



3(1 -x){l- vY 9(1 - x){l - yf 3(1 - a:)(l - yf 
304 40 704 224 

' 9(1 -x){l-y)^ 3(1 - a;) ^3(1 - yf{x + 1) ^ 3(1 - yf 

1760 560 2368 160 

3(l-y)4(x+l)+3(T^ + 9(l-y)3(x+l)"3(r:^ 
192 320 80 208 

+ 77 7777 TV - 7777 777 " 77 77 TT + 



(l-y)2(a; + l) 3(l-?/)2 {l-y)[x + l) 9(1 -y) 



9(a; + l)y ' ' ' ' V 9 'ix{l - yY 

32 4 4 26 80a: 

+ ^7 777 + 



3a;(l-|/)3 a;(l-|/)2 ^x{l - y) 9a; 3(1 - 1/)^ 
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+ 



32x2 64x2 
+ 
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3(1 3(l-y)2(a; + l) 3(1 - y)^ 
32 32 100 



+ 
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64 112 64 8 



+ 



3(1 3(l-2/)2(x+l) 3(1 -y)2 3(l-y)(x+l) 
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i/(0;x)/7(0;y) 



(100) 



sf°)(-p^-g^) = -{f + 
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4 



3(1 -y)3 9(1 -y)2 9(1 - y) 3(y + 1) 
16x 16x2 32x2 28x2 
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128 208 376 



3(1 -y)2(x + 1)3 3(1 -y)2(x + 1)2 9(1 _ y)2(a; + i) 
140 (31 128 
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+ 



3x2(1 -?/)3 27x2(1 -y)2 27x2(1 -y) 3x(l - y)^ 
352 1172 112 10 176x 

+ 7^^7 7^ - 77^-^ 777 - TT^^T T - 7^ + 
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20a;2 64 128 56 



+ 



9(1 -y) 9y(a; + l)4 9y{x + If 9y{x + lY 
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+ -— 


200 




28 

+ TTT 



64 



128 



+ 
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+ 



3(1 -y)2 9(1 -y) 3(y+l)2 3(y+l), 

16 32 28 20 



+ 



+ 



3a;2(l-y)4 3a;2(l - y)3 9a;2(l-y)2 9a;2(l - y) 
64 128 112 80 

+ 7r^^ TT - ITT-. + TTT-. TTT + 



3a;(l-7/)4 3a;(l-2/)3 ^x{l - 9x(l - y) 

64a; 128x 112a; 80a; 

+ 



3(l-2/)4 3(l-2/)3 9(l-2/)2 9(1-2/) 

16a;2 32a;2 28a;2 20a;2 

+ 



3(1 -y)4 3(1 -y)3 9(1 - y)^ 9(1 - y) 

64 128 112 
+ ^t:^ ^-T-. ^ + 



(l-x)(l-y)4 (l-x)(l-y)3 3(l-a;)(l-y)2 
80 64 64 



3(l-a:)(l-y) 3(1 - y)4(a; + 1) 3(1 -y)^ 
128 128 112 



3(l-y)3(a; + l) 3(1 - y)3 9(1 - y)2(a; + 1) 
112 80 80 

^(-1,0; x) 



9(1 -y)2 9(l-y)(a; + l) 9(1 - y) 

2y 



- - + 



+ 



2y 



3 3y(l-a;) 3y(a; + l) 3(1 - a;) 3(a; + 1) 

4 4 



3(1 -a;) (a; + l) 
20 4 



if(0,0,0;a;) 



+ 



-77 + 



3 3x2(1 - y)5 
2 2 



3a;2(l — y)4 a;2(l — y)3 3a;2(l — y)2 a;2(l — y) 
4 32 80 16 



3a;2(y + l) 3a;(l - y)^ 3a;(l - y)^ a;(l - y)3 
8 20 16 8 



3a;(l-y)2 3a;(l - y) 3a;(y + 1)3 x(y + l)2 

4 32,T 80a; 16a; 



3x(y + l) 3(1 -y)5 3(1 - y)^ (1 - y)3 

8a; 4a; 16a; 8,x 4a; 



3(1 -y)2 3(1 -y) 3(y + 1)3 (y + 1)2 3(y + 1) 

8a;2 20a;2 4a;2 2a;2 2a;2 

+ 



3(1 -y)5 3(1 -y)4 (1 - y)3 3(1 - y)2 3(1 - y) 
32 80 



4a;2 

+ 777 7T + 



3(y + l) (l-a;)(l-y)5 (l-x)(l-y)^ 



50 



+ 



48 8 8 

(1 -x){l- yf + (1 - x)(l - yf " (1 - x)(l - y) 
32 32 80 



+ 



3(l-y)5(x + l) ?,{l-yf 3(l-y)4(x + l) 
80 16 16 8 



3(1 (1-^)3(^ + 1) (i_^)3 3(l-y)2(x + l) 



+ 



16 32 

+ 



3(1 3(l-|/)(x+l) 3(1-1/) 3(1/ + 1)3 



16 



+ 



(t/+l)2 (I/+1) 

80 16 



^(0,0,0;y)+ - 



32 



3 3x2(1 -|/)5 

8 



+ 



3x2(1 -y)4 x2(l-y)3 3x2(1 -y)2 3^2(1 -y) 
128 320 64 32 



2>x{l — yY 3x(l — y)"^ a;(l — ?/)3 3x(l — y)2 
32 128,T 320x 64a; 32x 

T + TTT: 77-777: 77 + 77 77 4 



M^-y) 3(l-2/)5 ^{l-yY {l-y f 3(1 - 1/)2 

2 



32x 32x2 

+ 



80x2 iQ^2 

+ 77 77 + 



5X 



3(1-1/) 3(1 -y)5 3(l-|/)4 (l-y)3 3(l-|/)2 

2y 



8x2 2 

+ 



2 ^ 2y 



3(1 -y) 3?/(l-x) 3?/(x + l) 3(1 -x) 3(x + 1) 

128 320 192 

+ 



(1 - x)(l - yf (1 - x)(l - yY (1 - x)(l - yf 

32 32 4 

+ 



(l-x)(l-y)2 (l-a;)(l-y) 3(1 - x) 
128 128 320 

+ 7777 777 77 + 



3(l-y)5(x + l) ^{l-yf 3(1-|/)4(,T+1) 



320 



+ 



61 



+ 



64 



+ 



2>{l-yf (l-z/)3(x+l) {l-yf 3(1 - y)2(x+l) 
32 32 32 4 



3(1 -y)2 3(l-y)(x+l) 3(1 -t/) (x + 1) 

x//(0, 0; ,)//(0; x) + + ^^.^ 

4 2 2 32 



X 



+ 



X2(l-y)3 3x2(1 -7/)2 3x2(1-7/) 3x(l - y)^ 

80 16 8 8 

+ 7^^7 77 - -7Z 77 " T^^T 77 + 



?>x{l-yf x{l-yf 3x(l-?/)2 3x(l - y) 

32x 80x 16x 8x 

+ 77 77 + 



2>{l-yf 2>{l-yf {l-yf K^-yf 

•2 20x2 4x2 2x2 



5X 



5X 



3(1 -y) 2>{l-yf 2>{l-yf {l-yf 2>{l-yf 
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2x^ 



32 



^0 



3(1-1/) {l-x){l-yf {l-x){l-yY 
48 8 



+ 



{l-x){l-yf {l-x){l-yY {l-x){l-y) 

32 32 80 

+ TTT rr + 



3(l-y)5(x + l) 3(1 -y)5 3(l-y)4(x + l) 
80 16 16 



+ 



3(1 -y)4 [l-yf{x+l) [l-yf 3(l-|/)2(x+l) 



?>{l-yY ?>{l-y){x+l) ?>{l-y))[\^^ ^'^'^'"^^ 
-G(-l/|/, 0, 0; x)) + (G{-y, 0; a:) + G{-l/y, 0, a:))i/ (0; y) 

+ (G(-l/y; x) - G(-y; x)) (i/ (0, 0; y) + 3C(2)) 
+2//(0; x)//(l, 0; y) - 2H{-1, 0; x)//(0; y) 



-6i^(-l,0;y)i^(0;a;) 



;ioi) 



B Asymptotic Expansions 

In this Appendix, we provide the asymptotic expansions of the auxihary functions 
introduced in the paper. For the functions related to the self-energy and vertex 
correction, that depend on a single kinematic invariant, we give the expansion in 
the case in which the kinematic invariant is significantly larger than the electron 
mass m. For the functions B^^^^ and B^^^ {i = 1,2,3), related to the box graphs, 
that depend on two kinematic invariants, we discuss the expansion in two different 
cases: 

i) both arguments of the functions B^^^^ and B^'^''^ are larger than m^; no specific 
hierarchy between the two arguments is assumed; 

ii) both arguments are much larger than the electron mass, but the first (second) 
argument is also much larger than the second (first). 

Employing these results, in combination with the expression for the cross section 
in terms of the auxiliary functions, it is possible to obtain the Bhabha scattering 
cross section in the kinematic regimes in which s t u rn? , s t ^ u ^ m^, 
and s ~ ^ t ^ m^. 

The expansions of the auxiliary functions are truncated in such a way that the full 
cross section can be recovered up to terms of order /F? (where = —s, —t. —u) 
with » m\ In particular, for the functions n^^'), nf \ fI''\ f[''\ F^''^ and 
^2^''' one needs only the term of order {m^ / R^)^ . On the contrary, for the functions 
Bf^ and Bf^ (^ = 1, 2, 3) one needs only the terms proportional to {w? / R?) ^, 
since these functions appear in the cross section always multiplied by w? / s or w? /t. 
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In the following equations we use the definitions 



m 



—J, LQ2 = ln^^^j, 



(102) 



In the expansions we will keep and both in the Euclidean region; the 
continuation to timelike values is trivial; if, for instance, — > — (s + ie), then 



TTl 

Lp2 In ( — I + ZTT . 



(103) 



B.l One- and Two-Loop Auxiliary Functions in the Small 
Mass Limit P^^ > 

The first region of interest is the one in which P'^,Q'^ ^ w?, {x, y <g 1), relevant for 
the large angle Bhabha scattering. We find: 



n(i^,o)^_p2) 



n(iM)(_p2^ 



9 + 3^-^'' 



5 

18' 



1 

12' 



,(l/,0). 



-P2) 



1 +Lp2 , 

-1 + ic(2) - \lp. - \lU , 
0, 



(1^-1)/ 



(iz,o). 



p2 



= -16— Lp2-8-^L 



jp2 



jp2 



(104) 
(105) 

(106) 
(107) 

(108) 

(109) 



p2 



20C(2) + 2Lp2 - 21% - 4:Lp2Lq2 + 2Lq2 + 2LJ2] 
10(^(2) + 2Lp2 — L?p2 — 2Lp2Lq2 + 2Lq2 + -^'Qzj 



16C(2) + 4Lp2 - 4Lp2Lq2 + 2LL 



B. 



(11,0). 



-P\-Q') 



p2 Q2 pi 

■16 —Lp2 - 8 —Lp2 - 16 ^ Lp2 



(110) 



;iii) 



20C(2) - 2Lp2 - 2L^2 - 4:Lp2Lq2 + 2Lq2 + 2L^2 

10(^(2) + 2Lp2 — Lp2 — 2Lp2Lq2 + 2Lq2 + Lq2 

32C(2) - 8Lp2Lq2 + 4LJ2 , (112) 
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^(l/,-l)^_p2^_g2) 



p4 

8 „ „^ Lp2 , 



P 



B. 



(U,0), 



-P\-Q') = 4— Lp2+4^Lp2 + 



-\-1L'n2 



(113) 

16C(2) + 4Lp2 - ALp2Lq2 

(114) 



and for the two-loop functions: 



n, 



(2/,0), 



(115) 



^ _ 1^(2) + ^Lp2 + ic(2)Lp2 + + , (116) 

108 4^^ ^ 216 6^ ^ 72 ^ 36 ^ ^ ' 



,(2i,0). 



-P2) 



0, 



(117) 



p2 
p4 



9 



p(2^,0)^_p2^_g2^ 



p2 - 
2 



40 r 8^ ■ 

— Lp2 + -Lp2Lq2 

20 4 

— Lp2 + -Lp2LQ2 

i.ii,,H4c,.).:(i4 



20 r 4 

— Lp2 + -Lp2Lq2 



(118) 

242 13 2 
— Lp2--Lp. 



3L|,2 In (l + 1^ 



2 
"3 

34 



9 



Lp2Lq2 — 2Zyp2ZyQ2H -Zyp2 Li2 



16 



+8C(2)Lq2 + -L^ 



■-L^2--Lp2 LQ2ln( 1 + ^ 



20 

p2 



In ( 1 + ^ 1 + -Lj, 



( 

3 H'p^ 



| + fc(2)+2C(2)ln(l + |; 



130 



5 
9" 



Lp2 — — Z/p2 



+-L|,2 In I 1 + 



p2 



--Lp2--Lp2LQ2 In 1 + 



p2 



20 2 f 

-—Lp2Lq2 -Lp2LQ2 + -Lp2 Li2 ( --^ j +4C(2)Lq2 



+2Lc 



25 



;^^^ + 3 



L^2 In ( 1 + ^ 



2 / g^x 2 

-3i^Q^Li2(--p,l+-Li3 



2 

p2 
p2 



;^.>,..C(.,.n(:.|:)-|.. 



54 



+^C(2)ivp2 + In f 1 



1 
3 



p2 



- -^Lp2Lq2 



--Lp2Lq2 In M + — j + -Lp2 Li2 ( -j^ 



--Lp2Ll 



+ -Ll2\Ti[l + ^\+-L 



p2 



23 



^Lq2 + ^((2)Lq2 + ^Lj 

-i7Q2Ll2 I — — 



9 



p2 



(119) 



P 



>2 r- 



40 , 8 , , 
40 



r20 



-^-^p2 + -^Lp2Lq2 



m2Q2 
p2 r 
2 



g -^P2 + -^Lp2Lq2 



(120) 



'68 14,,^, / g2\ 272 , 13, 



2 

p2 



2 
"3 



+^Lp2 In ( 1 + ^ 



p2 



- 3 4 , f 

-Lp2 - -Lp2Lq2 In I 1 + 



2 
9' 



40 2 4 / g2\ 20 

+8C(2)Lq2 + ^LJ. + ^L|. In (l + + ^^1^ 



— -LQ2lA'y I — 



p2 



5 
9" 



^p2 — — /yp2 



- + -C(2)+2C(2)ln(l + ^j-- 

In (1 + - - f In (1 + %) 

20 2 2./ 

^Lp2Lq2 — Lp2Lq2 + -Lp2Li2 ( 1 + 2Lq2 

+4C(2)Lq2 + + In (1 + 1^) + 

2, / Q^V 2,. / Q' 



68 80 
27^Y 
224 



C(2)+4C(2)ln (^1 + + lc(2)Lp2 



27 



2 
3' 



Lp2 + ^Lp2 In ( 1 + 



p2 



4 ( 
--Lp2Lq2 In M + — 



- -^Lp2Lq2 

8 2 4 

-Lp2Lq2 + -Lq: 
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B?'-'\-P^-Q') 



B. 



(21,0) 



+ -Lp2Li2 



32 

y 

10 



C(2) V + ^L|. 



+ 3L13 



p2 



p4 



p2 - 



20 r 4 

— Lp2 + -Lp2LQ2 



(121) 
(122) 



10 r '^T T 

— Lp2 + -Lp2LQ2 



10^ 2 

— Lp2 + -Lp2LQ2 



p4 



-L^2 In ( 1 + ^ 1 + -Lp2Lq2 



+-Lp2Lq2 In M + — j + -Lp2L^2 - - 
2. . / 16,, 23, 



2,. / 

"3 



2 

3" 



alio - 



p2 

(123) 



B.2 One- and Two-Loop Auxiliary Functions in the Limit 
> > and > > 



In this limit, the Eqs. (ITintill8ll2()fT^ remain vahd, while Eqs. fjll9ll2lll23|l must 
be replaced, in the case in which ^ ^ iv?, by: 

[34 40,,^, 82^ 2,, ,^ 8^ ^ 

^ + Y'^'>2j - ^^P2 + -C(2)Lp2 - -Lp2Lq2 



{21,0), 



-P\ -Q' 



+- 



-Lp2Lq2 + -Lq2 +—Q[2)Lq2 + y^^q2 + 9 Q- 

■Lp2Lq2 



50 32.,^, 260^ 10^2 2^3 
^ + YC(2)-— Lp2--Lp.--Lp. 



9 



26 



19 



2Lp2Ln2 + —Lq2 + 8C(2)Lq2 + — Lq2 + -L 





' 85 


^2p2 


162 



9 

2 , ^ 11 
3C(2) + ^LP. 



9 



-/^p2 + -Lp2Lq2 



"27^^^ 



56 



115 
2592 
13 . 



1 1 S 

_((2) + —Lp2 
6^^ ^216 ^ 



36 



-L 



p2 



— Lp2Lq2 



216 



J^n2 



36 



1 49^,^, 323 ^ 1 2 lr3 
108 + ¥^^2^ "^^^^"18^^^ "9^^^ 



29 



19 



17 



+ ^^2] I , (124) 



p4 



68 80^,^, 

27 + 



224^ 4^,,^ 20^ ^ 

— Lp2 + -4(2jLp2 - —Lp2Lq2 



8^ ^2 4^ 32,, 23^2 lOrS ' 

--Lp2LQ2 + -i.Q2 + yC(2)i-Q2 + — LQ2 + yi.Q2 



+- 



32 50,, , 308 , 
+ — C(2)-^Lp2 



7 



27 



27 



■I p2 



2 

22 



52 



-Lp2Lq2 



-2Lp2Lq2 + + 8C(2)Lq2 + y-^Q2 + -Lq2 





" 89 




162 " 




-Hi 








27 




' 71 


^2p4 


1296 



;C(2) 



27 



-^P2 — g-^P2 + n-^P^Lg 



9 



11 

108 



L p2 



-Lp2 



-Lp2Lq2 



11 1 

-ivr)2 — 



108 



54 9 ^ 27 



2 
9" 



Lp2 — - Lp2 — - I/p2 ^ -^P2 -f'Q2 



26 



31 



-Lp2Lq2+3Lq2 +AC{2)Lq2+—Lq2 + -Lq2 



(125) 



4 



34 40,,^, 82^ 2.,^,^ 8^ ^ 
^ + y^*^^) ~ 27 3C(2ji^p2 - -Lp2Lq2 



16 



23 



Lp2Lq2 + -Lq2 + — C(2)Lq2 + — L 



18 



9 V 



+- 



P2 



'2 4 12 1 

2 - 2C(2) - -ivp2 - -Lp2 - -Lq2 - -Lq2 

2 2 2 ^ 2 



8l"3^^2^ + 27^^^"9^^^ + 9^^^^«^ 



57 



24 



12 
23 



12 



--L 



C(2) + Yg-^-P" ^ Q^P^ + Lp2Lq2 - -Lc 



(126) 



In the case in which ^ ^ m^, Eqs. (jll9ll21ll23|) must be replaced by: 



(21,0). 



p2 



50 32 



224 
~27 



4 
3' 



10 



9 



+ -Lp2LQ2 - —Lp2LQ2 - -Lp2Lq2 + —Lc 



20,, 19 . 

+ -C(2) V + -Lj: 



9 



73 



9 



34 22 



130 



- + -C(2) - —Lp2 + -C{2)Lp2 - -L 



27 9 
1 



27 



9 



p2 



20 



+ -Lp2Lq2 - —Lp2Lq2 - -Lp2Lq2 + 2Lq2 



10 



25 



+—C{2)Lo2 + —Li2 + -L 



p6 


' 85 




162 


27 


p4 


1 




108 



:c(2) 



Lo2 - -Lc 



11 
27 



L p2 



gLp2 + -Lp2Lq2 



"g-^p2 + -Lp2Lq2 —Lp2Lq2 — -Lp2Lq2 

"Y^^Q^ + yC(2)^Q2 + yi^ga + ^Lj: 



(127) 



p2 



50 32^, , 254^ 4,, 10 . 4 o 



2 

J 

3 

20 



— Lp2Lq2 

19 



-Lp2Lq2 + —Lq2 



14 



+-C(2)Lq. + -L^. + -L^. 



58 



"34 22,,^, 130, 2 
.27 + - 2f + 3' 

2 



3C(2)L, 

+yC(2)LQ. 



ip2 — gL%2 



g-Lp2LQ2 - -Lp2Lq2 + 2Lq2 



f2 + g^Q2 



p6 


"193 




162 


29, 
+ 27^ 


p4 


"137 




108 



27^ 



- + 



107^,^, 385 , 8^,^,, 1. 
-9-C(2) - ^i^p^ + 3C(2)LP2 + -L^. 



13 r ,28 
18 ^ 



10 
"3 



— -^p2ivQ; 



I 2 

-Lp2Lq2 



(128) 



p2 



10 r 2 

— Lp2 + -Lp2Lq2 



10 r 2 

— Lp2 + -Lp2LQ2 



3 - 2C(2) - -ivp2 - -L^a + -Lp2Lq2 + -L(; 



3V 

34 40^,^, 82, 4^,^,, 
^ + yC(2)-^Lp. + 3C(2)Lp. g 



— J-'p2 



+ -Lp2Lq2 - -Lp2Lq2 - -Lp2Lq2 + -L 



3 
14 

+^C(2)Lq 



23 2 ^ r 3 



(129) 
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